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ABSTRACT 

We show the effects of AGN-driven outflows on the ejection of heavy elements using our cos¬ 
mological simulations, where super-massive black holes originate from the first stars. In the 
most massive galaxy, we have identified two strong outflows unambiguously driven by AGN 
feedback. These outflows have a speed greater than ^ 8000 km s“^ near the AGN, and travel 
out to a half Mpc with ^ 3000 kms“^. These outflows remove the remaining gas 3 per 
cent of baryons) and significant amounts of metals 2 per cent of total produced metals) 
from the host galaxy, chemically enriching the circumgalactic medium (CGM) and the inter- 
galactic medium (IGM). 17.6 per cent of metals from this galaxy, and 18.4 per cent of total 
produced metals in the simulation, end up in the CGM and IGM, respectively. The metallici- 
ties of the CGM and IGM are higher with AGN feedback, while the mass-metallicity relation 
of galaxies is not affected very much. We also And ‘selective’ mass-loss where iron is more 
effectively ejected than oxygen because of the time-delay of Type la Supernovae. AGN-driven 
outflows play an essential role not only in quenching of star formation in massive galaxies to 
match with observed down-sizing phenomena, but also in a large-scale chemical enrichment 
in the Universe. Observational constraints of metallicities and elemental abundance ratios in 
outflows are important to test the modelling of AGN feedback in galaxy formation. 

Key words: black hole physics - galaxies: evolution - galaxies: formation - methods: nu¬ 
merical - galaxies: abundances 


1 INTRODUCTION 

The importance of feedback from active galactic nuclei (AGN) has 
been underscored by the discovery of the relationship between the 
mass o f the central black hole (BH) and the mass of the host galaxy 
bulge dMagorrian et alJfTggsl : iKormendv & Ho|[2013ll . suggesting 
co-evolution of BHs and their host galaxies. This has already been 
indicated from the simil ar shapes between th e observed cosmic star 
formation his tory (e.g., MadaU-gLaL I ll996h and the quasar space 
density (e.g., ISchmidt et alJll995h . AGN feedback has been im¬ 
plemented in cosmological simulations, which pr ovided an excel¬ 
lent agreement with the M agorrian relation (e.g.. lDi Matteo et al.l 
l2008l : ]^iacki et ai]|2014ll a nd a better reproduction of t he cosmic 
star formation rates (e.g.. iTavlor & Kobavashil l l20l4 hereafter 
iTKldll . lVogelsberger et al.l2014 ). The [a/Fe] problem in early-type 
galaxies also requires AGN feedback, which plays an essential role 
in quenching star formation i n massive galaxies where supernova 
(SN) f eedback is inefficient dTavlor & Kobavashil 1 20151 hereafter 
iTKlSll . 

For the quenching mechanism, SN-driven galactic winds have 
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been proposed dLarson|[l974l : lArimoto & Yoshiilll9^ . but are not 
efficient enough for massive galaxies in hydrodynamical simu¬ 
lation s with dark matter (e.g.. iKobavashil[20051 : iKobavashi et al.l 
l2007h . There is observational evidence for galactic winds both lo¬ 
cally and a t high redshifts, from low-mass star-forming galaxies 
(e.g., M82. lOhvama et al ] 2002h to AG N-hosting massive galax¬ 
ies (e.g., Centaurus A, Kraft et aT]|2009t) . Supernova-driven winds 
typically have velocities of a few 10^ — 10® km s~® and an out¬ 
flow rate comparable to th e star formation rate (e.g.. lHeckman et ^ 
120001 : jPettini et^ l2000h . while AGN-driven winds show much 
higher velocities and outflow rates. The presence of outflows in 
luminous q uasars and A GN has been evidenced in broad absorp¬ 
tion lines dUyndsiIl967h . The nearest quasar Mrk 231 shows a 
multi - phase outflow con taining ionized, neutral ( Ru£ke_&Willc^ 
20ljJ^ngrt^ 20141) . and molecular gas dFeruglio et al.ll201tt 


Cicone et alJl2012h . which may be explained by a bipolar outflow 


and an accretion disk. Winds are also seen in Seyfert galaxies (e .g., 

iTombesi et alj201lLl2013l: [Pounds & King||2013i:|PTOndsll2014l) . in 

which velocity probably dep ends on the distance from the central 
BHs on ~pc to ~kpc scales jTombesi et al.[[20r^ . and in other ul- 
traluminous infrared galaxies with velo cities depending n ot on star 
formation rate but on AGN luminosity l lSturm et al.[[201 ih . At high 
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redsh ifts, AGN-hosting gal a xies also show evidence of outflows 
(e.g., iNesvadba et alj 1200 Sl : iHarrison et alj I2OI2I : iMaiolino et all 
I2OI2L for z = 6). These suggest that AGN feedback should be ac¬ 
tivated at high redshifts. This is the case in our AGN m odel where 
super-massive BHs originate from the first stars ( ItK14|) . 

Galactic winds are the most pl ausible origin of mass- 
metallicity relations of galaxies (e.g., iTinslevI Il980ll : low-mass 
galaxie s have lower metallicit ies becau se of larger metal loss by 
winds jKobavashi et alj 120071 see also lOave et aD 1201 ih . These 
metals are distributed to the circumgalactic medium (CGM), intr¬ 
acluster medium (ICM), and intergalactic medium (IGM). The de¬ 
tailed mecha nism of this metal transport has not been well stud¬ 
ied yet (e.g., iBarai et ^ |20l4) . Direct o bservations of outflows 
from AGN (e.g., Kirkpatrick et al. I l201lh and of chemical abun¬ 
dances near AGN will provide another strin gent constraint on the 
modelling of AG N i n galax y simulations l l Germain et alj|200^ : 
iBarai et aObOl ih . In ItKISI we showed that the mass-metallicity 
relation is preserved when AGN feedback is included, even though 
the mass of metals in massive galaxies is reduced; in this Letter we 
will investigate how the matter and heavy elements are redistributed 
from galaxies to the CGM and IGM by AGN-driven winds. 

In this Letter, we provide a new method to identify the out¬ 
flowing gas in chemodynamical simulations (Section |^. We ap¬ 
ply this to the most massive galaxy in our cosmological simula¬ 
tions, and investigate two unambiguously AGN-driven outflows. 
We then describe their impact on the metallicity of the host galaxy, 
CGM, and IGM in Section [S] Finally, we present our conclusions 
in Section [4] The resulting spacial distributions of metallicity, and 
the IGM metallicity, will be shown in a future work (Taylor & 
Kobayashi 2015, in prep.). 


2 DEFINITION OF OUTFLOWS 

Our simulations use the GADGET-3 code dSpringell l2005h . 
augmented with physical processes relevant to galaxy fo rma- 
tion and evolution: star formation jKobavashi et al.l l2007h . en- 
ergy feedback and chemical enrichment from T ype II and la SNe 
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) and hvpernovae 

Kobayashi et alj|2006l; Kobayashi & Nakasatol 

201 ib; BH forma- 


tion, growth through gas accretion and mergers, and energy feed¬ 
back jTK14ll : heating from a uniform, evolving UV background; 
radiative gas cooling. The simulation box is 25h~^ Mpc on a side, 
with initial conditions for 2 x 2 40^ pa rticles that produce a central 
cluster by the present day (see jTKlSh for more details). We use a 
gas gravitational softening length e = 1.125h~^ kpc. 

We detect inflows and outflows by tracing gas particles for the 
lif etime o f the simulated galaxy. In this Letter, we focus on galaxy 
A jTKlSl) , which is the largest galaxy at z = 0 in our cosmolog¬ 
ical simulations. This galaxy sits at the centre of the cluster, and 
grows to M* = 6.8 x IO^^Mq, Mtot = 9.5 x 10 ^^ Mq, and an 
effective radius of kpc by 2 = 0 (see Table 1 of lTK15l for 

Other quantities). We first look at the net radial galactocentric mo¬ 
mentum (measured in the galaxy’s rest frame) of gas particles to¬ 
wards the outskirts of the galaxy, specifically at galactocentric radii 
20h~^ kpc < r < 25/z“^ kpc, for the lifetime of the galaxy. This 
range was chosen to be close enough to the galaxy that gas particles 
passing through the shell could be associated with the galaxy, but 
sufficiently distant that the orbits of gas within the galaxy would 
not affect it. This quantity is shown as a function of time for both 
simulations with and without AGN feedback in Fig. \I\ Negative 
values indicate a net inflow of gas; at late times, the galaxy without 



Figure 1. Net radial momentum of all gas particles in a shell 20 — 2bh~^ 
kpc from the centre of the largest galaxy (at 2 : = 0) in the simulation across 
cosmic time. Positive and negative values indicate a net outflow and inflow 
of gas through the shell, respectively. This value is shown for the simulation 
with (without) AGN feedback in the top (bottom) panel. In both panels, blue 
vertical lines indicate the occurrence of major mergers. The blue arrows 
indicate the times that the large outflows are first identified by the method 
of Section [ 2 ] 



Figure 2. Radial distribution of the net radial momentum of outflowing gas 
particles (see text for definition), within 500/i~^ kpc of the most massive 
galaxy in the simulation with AGN. We show this distribution at several 
different times (labelled relative to 9.201 Gyr after the Big Bang), allowing 
the progression of the outflow away from the galaxy to be seen. 
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Table 1. Total and metal masses, and abundances, of the outflows, ISM, CGM, and IGM. 


Region 

z 

Mg/10i“MQ 

Mz/1O®M0 

[Fe/H] 

[0/H] 

[O/Fe] 

Outflow 

0.41 

12.3 

268.5 

-0.78 

-0.96 

-0.09 

Outflow 

0.31 

8.6 

196.6 

-0.75 

-0.94 

-0.18 

ISM 

0.41 

0.07 

35.9 

0.62 

0.39 

-0.22 

ISM 

0.31 

0.08 

21.7 

0.17 

0.17 

0.14 

CGM (AGN) 

0.00 

87.3 

2302 

-0.67 

-0.78 

-0.14 

CGM (No AGN) 

0.00 

216 

4292 

-0.73 

-0.96 

-0.20 

IGM (AGN) 

0.00 

2.08e4 

1.14e5 

-1.77 

-1.72 

0.00 

IGM (No AGN) 

0.00 

2.02e4 

8.62e4 

-2.12 

-1.92 

0.16 


AGN (bottom panel) shows prolonged accretion of gas, punctuated 
by short-lived, SN-driven winds. The same is not true of the galaxy 
when AGN feedback is included (top panel); at late times, there are 
a number of long-lived events in which the galaxy loses gas, most 
clearly between 10.5 and 12 Gyr. In both panels, blue lines indicate 
the occurrence of major mergers (defined as having a mass ratio of 
1:3 or greater). 

Although disentangling the contributions from AGN and stel¬ 
lar processes is not trivial, two very massive outflows at f 9.2 
Gyr (z = 0.41) and t ^ 10.1 Gyr (z = 0.31) are unambiguously 
AGN-driven, and so we further analyse these. In order to define the 
outflows, we determine which gas particles form them. Since we 
are searching for AGN-driven outflows, we expect their constituent 
particles to be travelling radially from the galaxy centre. Therefore, 
we select those particles within 500/i“^ kpc of the galaxy whose 
radial momentum exceeds 95 per cent of its total momentum, p, i.e. 
p ■ r > 0.95p. Choosing a slightly higher or lower threshold does 
not affect our results in the following sections. Fig. [Ogives a visu¬ 
alization of the outflow at f ~ 9.2 Gyr from the centre of galaxy A, 
showing the radial distribution, at several epochs, of the particles 
obtained by the above method. Times are given relative to the time 
of the snapshot immediately preceding the outflow’s identification, 
9.201 Gyr after the Big Bang (black line). The progression of the 
outflow from the galaxy into the IGM is clear, and it can also be 
seen to become broader and more massive as it sweeps up other 
gas particles. 

The total masses and metal masses of the AGN-driven out¬ 
flows when they are first identified (red line of Fig[3 are given 
in Table [T] In order to decouple the outflow caused by the AGN 
and the (AGN-enhanced) SN wind, we ignore the particles iden¬ 
tified from the snapshot immediately preceding it (i.e. the black 
line in Fig. |2j. Those particles that remain are defined to be the 
AGN-driven outflow; they are hot (T > 10^ K) and low-density 
(~ 10®mn m“®). On scales > 100 kpc, they expand i sotropically 
with a bubble-like structure, different to the findings of ICosta et al.l 
( l2014h who use the AREPO code; we leave a more detailed investi¬ 
gation of the outflow structure to a future work. Compared with the 
total gas mass 4.4 x IO^Mq in R200, the outflow mass is rather 
small. However, over the nex t several 100 Myr they travel further 
into the IGM, heating ( lTK15h and chemically enriching it (Section 
0 . 

We also estimate the velocity of these outflows. To each of the 
radial distrihutions in Fig.l^we fit a Gaussian function in order to 
estimate its center, from which we calculate an average speed at 
various stages of the outflow. These velocities are 3900 and 3400 
kms“® near the galaxy for the outflows at z = 0.41 and 0.31, 
respectively, falling to ~ 1000 kms“® at larger scales and later 
times. With values for the mass and speed of the outflows, we can 
also estimate their kinetic energies (albeit with large uncertainties) 
to be ~ ]^q 60-61 gj.g Note that the physical scale of our wind is 


much larger than for the available observations (Section[TJ, and our 
estimates of velocity and energy are larger than the observational 
estimates. 

With our usual timestep between simulations snapshots (~ 
25 — 45 Myr at these redshifts), these outflows are identified when 
they are already far from galactic centre (> 100 kpc). Therefore, 
we re-run the same simulation from the snapshot immediately pre¬ 
ceding the outflow at 9.201 Gyr (i.e., from the black line to the 
red line in Fig. 0, producing a larger number of snapshots with 
timesteps of ~ 175, 000 yr. This allows us to trace the initial phase 
of the outflow, which is necessary to understand its origin. The 
super-massive BH in this galaxy undergoes a period of accretion 
at 75 per cent its Eddington rate, which heats a small amount of 
gas close to the BH. On such a small scale, the gas is not dis¬ 
tributed homogeneously and isotropically around the BH, and so 
it is not straightforward to estimate the velocity. If we measure the 
initial velocity with a kernel-weighted average of the radial veloc¬ 
ity of the 72 gas particles nearest to the BH, we obtain 8000 
kms“^, which is comparab le to observations of AGN - driven out¬ 
flows on very small scales jTonflresi_e^J l20l^ . 1201 ll ; iKing et all 
I 2 OI 3 I ; iTombesi et '^120131 ; IPoundsI 20141) . We can also estimate 
the mass outflow rate to be ~ 400Mq yr~^, which lies within 
the observed range IQMp, vr~ ^ < Mo ut < lO®M 0 yr ® (e.g., 
iFeruglio et al.l201(il ; ISchonell et alj20l^ . The outflowing gas soon 
slows and entrains more gas, eventually stalling ~ 100 kpc from 
the BH by f = 9.245 Gyr, when the re-run simulation ends. The 
outflow seen on larger scales in the next snapshot then suggests that 
there are at least 2 separate accretion events, separated by at most 
88 Myr, which power the outflow seen at « = 0.41 on the largest 
scales. The AGN energy produced by this accretion can account for 
the kinetic energy of the large-scale outflow. 


3 METAL CONTENT OE OUTELOWS 

Our simulations also track the elemental abundances of gas parti¬ 
cles, allowing us to see how the outflows affect the chemical en¬ 
richment of the host galaxy, the CGM, and the IGM. We show in 
Fig.0 the net flux of metal mass through a shell at 20h ^ kpc 
< r < 25/i“® kpc from the centre of the galaxy, analogous to the 
radial momentum flux of Fig. [T] As in Section^ there is a trend for 
net accretion of metals when AGN feedback is not included, but a 
net loss of metals with AGN. 

This can be seen even more clearly in the top panel of Fig. [4] 
in which we show the metal mass in the central comoving 
kpc of the galaxy, normalized hy the total metal produced by stars 
(by 2 = 0), as a function of time. With AGN, the galaxy has 3.4 x 
10 ®®Mq of stars in its central kpc, while without AGN the 

stellar mass is 1.2 x 10^^ M^ . Altho ugh the star formation rate 
is lower at all times (Fig. 4 of lTK15h . initially the metal mass is 
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Figure 3. Net flux of metal mass of all gas particles in a shell 20 — 25h”^ 
kpc from the centre of the largest galaxy at 2 : = 0 in the simulation, across 
cosmic time. Positive and negative values indicate a net outflow and inflow 
of metals through the shell, respectively. This value is show for the simu¬ 
lation with (without) AGN feedback in the top (bottom) panel. The blue 
arrows indicate the times that the large outflows are first identified by the 
method of Section[2| 



Time [Gyr] 

Figure 4. Evolution over cosmic time of the ratio of the metal mass to the 
total produced metals (= yM*, with the net stellar yield y = 0.02) in the 
central comoving kpc of the galaxy. The black (red) line indicates 

the simulation with (without) AGN feedback. The blue an'ows indicate the 
times that the large outflows are first identified by fhe method of Section[2] 


larger with AGN feedback, since star formation and SN feedback 
is suppressed early in the life of the galaxy, and therefore a lower 
fraction of metals are ejected compared to the simulation without 
AGN feedback in which SNe alone are able to eject some metals. 
At later times, after the BH has grown to be super-massive, strong 
AGN feedback removes a much larger fraction than SNe alone. The 
blue arrows in Fig. |4] denote the times that the large AGN-driven 
outflows occur, as determined by the method of Section|2 It can be 
seen that they correspond exactly with sudden drops in the metal 
content of the galaxy. In fact the values of these drops match very 
closely with the values given in Table[T] obtained by identifying the 
outflowing particles at much greater radii, lending credence to the 
method in Section 

Similar drops are seen for the gas mass and the gas fraction. 
At ~ 8 Gyr, this galaxy contains a significant amount of gas (a 
few x10®Mq in 10h~^ kpc). The gas mass suddenly decreases 
due to the outflow at 9.2 Gyr, but recovers quickly in only a few 
X100 Myr, which suggests that cold streams are not disrupted for 
long, even with this powerful outflow. However, by tracing the ID 
numbers of gas particles, we find that most of the outflowing gas at 
10.1 Gyr does not fall back to the centre later. At > 10 Gyr, the gas 
mass becomes lower than 10® M©, which results in a gas fraction 
below one per cent, significantly lower than five per cent in the case 
without AGN feedback. 

The metallicities of the outflowing gas are much lower than 
the average gas-phase metallicity of the galaxy in 10h~^ kpc (inter¬ 
stellar medium (ISM) in Table[TJ. This is because the gas expelled 
by the BH originates outside the galaxy from cold flows that have 
undergone much less chemical enrichment; by tracing the particles 
back through the simulation to high redshift, we find that ~ 97 per 
cent of the gas that constitutes the z = 0.41 outflow originates in 
cold flows. The ISM metallicity increases slightly from 0.33 to 0.37 
at the « = 0.41 outflow, due to the removal of this low-metallicity 
gas. 

Some of the ejected material accretes back onto the galaxy, 
and is present in the second large outflow 860 Myr later. At 
f > 12 Gyr, the metal mass of the ISM increases (Fig. 0. We can 
identify this origin by tracing the particles, and find that it is due 
to both in-situ star formation within the galaxy and mergers with 
gas-rich satellites. As a result, the present gas-phase metal fraction 
is only a few per cent lower than the case without AGN feedback; 
this is the reason why AG N feedb ack does not affect very much the 
mass-metallicity relations ( iTKlSh . 

The AGN-driven outflows remove a significant fraction of the 
galaxy’s metal mass, transporting it to the CGM and the IGM. 
Table [T] gives the present metallicities of the CGM, measured in 
r < 150h~^ kpc, and of the IGM, which i s calculated for all gas 
particles that do not belong to galaxies (see lKobavashi et ^120071 
for the definition). In this simulated galaxy, the oxygen mass of 
the CGM is 8.1 x 10^ with AGN, and 1.4 x 10 ®Mq without; 
iTumlinson et alj | |201 ih estimate a lower limit of ~ IO'^Mq from 
observations. However, the fraction of CGM metals relative to the 
total produced in the Galaxy A is higher with AGN feedback (17.6 
per cent and 13.2 per cent with and without AGN feedback, respec¬ 
tively), which results in the higher metallicity, log ZjZ^ = —0.78. 
With AGN feedback, 18.4 per cent of the total produced metals in 
the simulation box is in the IGM at z = 0, compared with only 9.6 
per cent when it is not included. The IGM metallicity is also higher 
than the case witho ut AGN, and log Z/Z p, = —1.78; this is higher 
than the estimate bv ISchave et alj ( 2003h . but there is a large varia¬ 
tion in metallicities wi thin the simulation box . Note that, different 
from the simulation in iKobavashi et al] ( l2007h . our simulation has 
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a central concentration of galaxies, which may cause the relatively 
high metallicity of the IGM. 

Finally, we find a difference in the elemental abundance ratio 
in the lower redshift outflow (see Table[T}; [0/Fe] ratio of the out¬ 
flowing gas is —0.18, which is lower than for the ISM in the galaxy 
([O/Fe] = 0.14). The opposite trend is seen in the higher redshift 
o utflow. This may c orrespond to the selective mass-loss discussed 
in iKobavashil (l2005h . Even after star formation is quenched by the 
first outflow, Fe is produced by SNe la with a long delay. The sec¬ 
ond outflow may be collecting this SNe la contribution. On the 
other hand, S N-driven outflows s hould contain higher [a/Fe] ra¬ 
tios as seen in iMartin et alj (l2002h . Observational estimates of the 
metallicities and elemental abundance ratios of outflows are very 
important to put further constraints on the modelling of AGN feed¬ 
back. 


4 CONCLUSIONS 

We have shown that AGN feedback drives massive outflows from 
the most massive galaxy in our cosmological simulations. These 
outflows, together with metal-rich SN-driven winds, remove gas 
and metals from the hos t galaxy, causing the CGM and the IGM 
to become hotter dTKlSh and display greater chemical enrichment 
compared to the simulation without AGN. Nonetheless, the mass- 
metallicity relation of galaxies is not affected very much by AGN 
feedback. 

With AGN feedback, it is clear from Figs. [T] and that the 
net rate of inflow of both gas and metals into the galaxy is reduced 
across cosmic time. Focusing on the two most massive outflows 
detected by the net flows, we visualize the motion of outflowing 
gas (Fig.O and quantify the AGN-driven mass and metal ejection 
(Table 1). By tracing the particles, we find that most of outflowing 
gas originates in metal-poor cold flows, and that most of metals 
in outflowing gas does not fall back to the centre of the galaxy 
by the present. We also find evidence for ‘selective’ mass loss of 
iron relative to oxygen (Section]^. AGN feedback is an important 
process not only to quench star formation at late times in massive 
galaxies, but also to explain chemical enrichment in circumgalactic 
and intergalactic medium. 
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